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Abstract 

A tight binding model is introduced to describe the 
strong interaction limit of excitonic ordering. At sto- 
ichiometry, the model reduces in the strong coupling 
limit to a pseudo-spin model with approximate U(4) 
symmetry. Excitonic order appears in the pseudo- 
spin model as in-plane pseudo-magnetism. The U(4) 
symmetry unifies all possible singlet and triplet order 
parameters describing such states. Super-exchange, 
Hunds-rule coupling, and other perturbations act as 
anisotropics splitting the U(4) manifold, ultimately 
stabilizing a paramagnetic triplet state. The ten- 
dency to ferromagnetism with doping (observed ex- 
perimentally in the hexaborides) is explained as a 
spin-flop transition to a different orientation of the 
U(4) order parameter. The physical mechanism fa- 
voring such a reorientation is the enhanced coher- 
ence (and hence lower kinetic energy) of the doped 
electrons in a ferromagnetic background relative to 
the paramagnet. A discussion of the physical mean- 
ing of various excitonic states and their experimental 
consequences is also provided. 



1 Introduction 

The unexpected discovery of high-T c itinerant fer- 
romagnetism in doped hexaborides has re-ignited 
interest in the problem of excitonic ordering near the 
semiconductor-metal transition.^, ||| Excitonically 
ordered states are characterized by an off-diagonal 
order parameter describing pairing between conduc- 
tion electrons and valence holes. Early theoretical 
work by Volkov et. aJ.Q anticipated the emergence 
of ferromagnetism on doping such an excitonic state. 
These authors considered the limit of nearly nested 
overlapping conduction and valence bands with weak 
repulsive electron-electron interactions. In this limit, 



the problem can be approximately cast into a form 
nearly identical to BCS theory, and studied using 
the techniques of mean-field theory. Although this 
work (and some subsequent recent studies ||) suffers 
from the important physical mistake of neglecting 
the instability to phase separation, ferromagnetism 
remains nevertheless a generic feature in a corrected 
treatment. (6|, ^] 

While the appearance of ferromagnetism in the 
weak-coupling limit is encouraging, it is far from 
a conclusive and complete theoretical explanation 
for the experiments. First, Coulomb interactions 
in the hexaborides are not particularly weak, and 
most likely are comparable to the Fermi energy 
and band overlap. Second, the above explanation 
appears to hinge on the first-order nature of the 
excitonic to normal (E-N) transition in the BCS 
limit. While this feature, mathematically analogous 
to the first-order transition to the normal state due 
to pair-breaking by an external Zeeman field in a 
superconductor, [||, ||] is present in the nested mean- 
field limit, there do not appear to be any general the- 
oretical grounds mandating this behavior more gen- 
erally. Moreover, the universality of the experimen- 
tal results, now observed in a large number of differ- 
ent compounds (Cai_ x La x Bg, BaBg, Cai-^Ce^Bg, 
SrBg, . . . )||ro||, argues for the robustness of the phe- 
nomenon. 

To determine whether excitonic ferromagnetism is 
indeed more general than its weak-coupling theoret- 
ical basis, we consider here the completely oppo- 
site strong-coupling regime. This is not expected 
to be directly applicable to the hexaborides, as 
these materials are most likely best described by 
an intermediate-coupling model. Nevertheless, many 
useful insights are gained from this complementary 
limit. As usual, the principle assumption of the 
strong-coupling limit is the dominance of potential 
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over kinetic energy. This is achieved concretely us- 
ing a tight-binding model (see Eqs. |I||7|, in Sec. §), 
in which the conduction and valence bands of the 
conventional continuum theories are replaced by lo- 
calized a and b orbitals, respectively. The analog of 
band gap in the continuum model is the level split- 
ting Eq = E a — Efr > 0. The order parameter char- 
acterizing excitonic ordering is then a matrix in spin 
space: 



A 



a/3 



(1) 



where a a creates an electron with spin a =|, j in the 
a orbital, and bp annihilates an electron with spin (3 
in the b orbital. Excitonically ordered states thus 
have some partial occupation of the nominally ex- 
cited a states, result of Coulombic repulsion. In 
general, A Q/ g is a proper order parameter (i.e. one 
which characterizes a spontaneously broken symme- 
try) if the a and b orbitals have different symmetries. 
In this paper, we consider a "minimal model" with 
this property, comprised of one a and one b orbital 
per unit cell - see Fig. [l|. This mimics the situation 
in the hexaborides, for which the conduction and 
valence states also transform as different representa- 



tions of the cubic point group[ll]. Because of com- 
plications arising from orbital degeneracy, however, 
the appropriate representations for the hexaborides 
are three-dimensional rather than scalar. We defer 
the possible complications arising from these addi- 
tional degrees of freedom to a future investigation. 

As for the more familiar Hubbard model (see, e.g. 
Ref. [ 12 1 ) , the problem simplifies somewhat in the 



strong-coupling limit. Considering first the undoped 
system (half-filled = two electrons per unit cell), we 
obtain a novel quantum pseudo-spin model (Eqs. ||- 
p^ , Sec. ||). Within this model, the excitonic insula- 
tor (EI) appears as an intermediate state separating 
not a metal and a semiconductor but a Mott insu- 
lator and a semiconductor (or band insulator). In 
some respects, the behavior is argued to be quite 
similar to that of a quantum spin-1/2 XXZ antifer- 
romagnet in a magnetic field, with excitonic ordering 
analogous to XY antiferromagnetism. The "spins" of 
the model, however, can take on five distinct states 
per site: one singlet state with both electrons in the 
lower-energy b orbital, and four different spin states 
with one a and one b electron. This is in contrast to 
the two states of a single spin-1/2 particle. 

In the strong coupling limit, this large Hilbert 




Figure 1: Imaginative illustration of a model for 
which the tight-binding description employed phe- 
nomenologically here directly applies. Red circles 
and blue-l- green crosses represent s and d xy orbitals, 
respectively. 

space is "unified" by several approximate symme- 
tries valid at different energy scales. At the largest 
energy scales this is an enormous SU(4) group, cor- 
responding to arbitrary complex rotations of the 
four components of Aq,^. The approximate SU(4) 
symmetry fully unifies all possible excitonic states, 
including singlet, triplet, and singlet-triplet coexis- 
tences. These are described by the general decom- 
position 



A 
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A s l + A t ■ <? 



(2) 



where A s , At are the singlet and triplet order param- 
eters, and 2 and <r are the 2x2 unit and Pauli ma- 
trices in spin space, respectively. A system with ap- 
proximate SU(4) invariance contains the germ of fer- 
romagnetism, since several possible excitonic states 
(those with non-zero ReA s Ajf and/or ImAj A A|) 
give rise to net exchange fields, and hence a mag- 
netic moment. SU(4) symmetry implies that these 
states are low in energy. At intermediate energies 
the SU(4) symmetry reduces to an SU(2)xSU(2) in- 
variance, which reflects separate spin rotations of 
the a and b electrons. The latter is a symmetry 
of the conventional continuum models of Els, and 
transforms the order parameter in a "chiral" man- 
ner: A ulAU R , where U L and U R are SU(2) 
matrices. Finally, further weak interactions reduce 
this to a simple SU(2)xZ2 symmetry at the (very) 
lowest energies. 
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These symmetry considerations underly the sim- 
ple physical mechanism for ferromagnetism eluci- 
dated here ||l3|| . The dominant tendency imposed 
by Coulomb interactions is to excitonic ordering. 
With approximate SU(4) symmetry, however, the 
"orientation" (form of A a p) of the order parameter 
is nearly free and fixed only by weak "anisotropy" 
terms. In the undoped material, these anisotropies 
favor a simple paramagnetic triplet state. Doping 
introduces additional exchange energy contributions 
that modify the anisotropy, causing A a p to "flop" 
into a different orientation with a ferromagnetic mo- 
ment. In the present model, the excitonic order in 
the ferromagnet is of non-collinear triplet type, in 
which 

A s = 0, A t A A* + 0. (3) 

As shown in Sec. ||, in addition to ferromagnetic 
magnetization, this state has additional spatially- 
varying local static moments and spin currents trans- 
verse to the axis of net magnetization. The transi- 
tion to this state from the paramagnet is generally 
first order, and therefore coincides with a jump in 
the electronic density. Since experiments are per- 
formed at fixed charge density (dictated by the con- 
centration of dopant ions), the intermediate "forbid- 
den" range of dopings can be accommodated only by 
phase separation. With long-range Coulomb interac- 
tions included, macroscopic phase separation is im- 
possible, and charge domain formation is expected, 
as pointed out already in Refs. || [7j . 

The detailed demonstration of this behavior with 
doping is non-trivial. As for many other strongly 
correlated systems, the problem of doping is much 
more difficult than that of the stoichiometric Mott 
insulator. Indeed, as the EI state lies intermediate 
between band and Mott insulators, doping the EI is 
a sort of interpolation between doping a conventional 
band insulator and doping an antiferromagnetic in- 
sulator. The latter problem is of course at the crux 
of the physics of high-temperature superconductiv- 
ity, so that perhaps the experimental and theoreti- 
cal insights gained in the hexaborides will be helpful 
more generally. At any rate, doping the EI can be 
shown by very simple arguments to favor ferromag- 
netism in strong coupling. Essentially, the physics of 
this behavior is similar to the "Nagaoka effect" (13] 
in a doped antiferromagnet - ferromagnetic align- 
ment of the excitonic order parameters allows for 
more coherent propagation of the doped electrons, 



and hence a lowering of their kinetic energy. This 
mechanism is actually stronger in the EI than in the 
antiferromagnet, because of the global coherence of 
the excitonic condensate, and the near degeneracy 
(due to approximate SU(4) symmetry) of ferromag- 
netic and paramagnetic states. To provide a con- 




Figure 2: Partial phase diagram of the strong- 
coupling model as a function of E G (half the bare 
splitting between a and b orbitals), and [/,, the chem- 
ical potential. The shaded region is not analyzed in 
this paper. Thick lines indicate the boundary be- 
tween the undoped region below (with two electrons 
per unit cell) and the doped region above (with more 
than two electrons per unit cell). The AF 2 and AF 3 
phases are antiferromagnetic Mott insulators with 
two and three electrons per unit cell, respectively. 
The BI state is the band insulator. Intermediate be- 
tween the BI and AF 2 phases are the excitonic insu- 
lator (EI) and an insulator with coexisting excitonic 
and Neel order (EI/AF). The FM, FM*, PFM, and 
PFM* phases are all ferromagnetic metals (see Ta- 
ble H for the differences between these states), while 
PM indicates a paramagnetic metallic phase. All the 
metallic states above exhibit excitonic order. 

crete demonstration of these ideas, the strong cou- 
pling zero temperature phase diagram of the model is 
calculated in this paper using a "free Fermi gas" ap- 
proximation. This approximation captures the most 
important single quasiparticle physics of electronic 
propagation in an excitonically-ordered background, 
but neglects interactions between these quasiparti- 
cles. For simplicity, we also assume a fixed ampli- 
tude, Tr At A = Aq/2, of the excitonic order param- 
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eter. The latter assumption is valid for weak doping, 
x <C 1, in which the orientation of the ordering is of 
paramount importance. Putting together the results 
of this calculation and the stoichiometric behavior, 
we arrive at the partial phase diagram in Fig.g. This 
is in agreement with the general expectations stated 
above. It should be stressed, however, that this anal- 
ysis of doping is far from exhaustive. More detailed 
investigations of both the weak and strong coupling 
limits are currently underway) 15 1 

The remainder of the paper is structured as fol- 
lows. In Sec. ||, we present a detailed exposition 
of the (simplest) tight-binding model capable of de- 
scribing excitonically-ordered states, and consider 
the limit of infinite interaction strength. The bulk 
of the paper is contained in Sec. ||, where the model 
is analyzed for large but finite interactions, focus- 
ing on the stoichiometric situation with two valence 
electrons per unit cell. For this electron density the 
model is insulating, but can sustain excitonic and 
other types of ordering. The properties of the model 
with doping are discussed in Sec. We conclude 
in Sec. || with a clarifying discussion delineating the 
physical properties of various possible excitonic in- 
sulators, and the relation of the results of this paper 
to the hexaborides. 

2 Tight-Binding Model 

2.1 "On-site" terms 

We consider a minimal model capable of exhibit- 
ing excitonic order, which contains two orbitals per 
unit cell, so as to give rise to two bands in a non- 
interacting limit (the actual situation in the hexa- 
borides is more complex, with orbital degeneracy 
leading to multiple electron and hole pockets). A 
strong-coupling limit is obtained by first considering 
only local interactions within a unit cell: 

H = J2Ea(4a i -blb l )-^[a\a t +blb l 

i 

where a, b are electron annihilation operators for 
the "conduction" and "valence" states, respectively, 
obeying {a ia ,a} jj3 } = {b ia ,tf jf3 } = SijS a p. Here and 
throughout the paper, we use Latin indices i, j, . . . to 
denote the lattice site, and Greek indices a, (3, ■ ■ ■ =t 
, I to denote the spin state. Labels will be sup- 
pressed and implicit wherever clarity allows. The 



parameters E G ,fj,,U,V describe the "band gap" (or- 
bital energy difference), chemical potential, on-site 
"Hubbard" repulsion, and nearest neighbor repul- 
sion, respectively, within the unit cell. 

A crucial feature of Hq is the absence of direct 
hopping between the a and b orbitals within the unit 
cell. For excitonic ordering to be well-defined, it is 
necessary at a minimum that the a and b states be 
distinguished by a discrete symmetry operation, e.g. 
parity. When this is the case, direct hopping be- 
tween these orbitals is prohibited. It may be helpful 
to imagine an artificial situation in which the a and 
b orbitals represent s and d xy orbitals on a single site 
of a square lattice (see Fig. |]). 

In this situation, a and b orbitals are orthogo- 
nal both on the same site and on nearest neighbor 
sites. An overlap is possible, though for next-nearest 
neighbor pairs, i.e. on a diagonal. In general, an ex- 
change interaction is allowed by symmetry, and takes 
the form 



H\ — —Jh ^2 Sia ■ Sib, 



(5) 



where 5,;, 



it- a 
-ajcra,, bib 



-b]ab 4 . Here and in 



the following, the Pauli matrices a act in the spin 
space. On physical grounds, a ferromagnetic ex- 
change (Jh > 0) is most appropriate due to Hund's 
rule effects for orthogonal orbitals. For pedagogical 
purposes, we may wish to consider instead the oppo- 
site antiferromagnetic sign for this exchange. From 
the discussion in Sec. 1, it is clear that an essential 
ingredient for excitonic ferromagnetism is the near- 
degeneracy of singlet and triplet states. To build this 
into the strong coupling model thus requires small 
Jh- For the majority of the paper, therefore, we will 
neglect Jh or treat it as a small perturbation. 

2.2 Infinite interaction limit 

The analysis of the strong-coupling limit begins by 
first considering the on-site Hamiltonian, H s [ te = 
Hq + H\, in the absence of electron hopping be- 
tween adjacent unit cells. This may be thought of 
as the analog of the U = oo analysis of an ordi- 
nary Hubbard model. In this case, the occupation 
of each orbital is a good quantum number, and the 
states can be straightforwardly enumerated. Assum- 
ing E G > V > 0, and at first also J a b = 0, one ob- 
tains the phase diagram shown in Fig. ^. For the 
present study, we are particularly interested in den- 
sities near two electrons per unit cell. Note that the 
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Figure 3: Strong coupling (ultra-local) phase dia- 
gram in the E G -\i plane, neglecting exchange and 
all inter-cell hopping and interactions. Regions with 
zero and four electrons per unit cell are not shown. 
In each phase, the lowest-energy states are pictured, 
with a orbital above and b orbital below. 

doping behavior (i.e. on increasing fj,) in this regime 
depends crucially on the relative strength of E G and 
U. In particular, for 2E G > U — V , the preferred 
charge Q = 2e state is one with both electrons in 
the lower orbital, corresponding to the band insula- 
tor. For U — V > 2E G , by contrast, the two-electron 
ground state has one electron in each orbital, and 
hence a net spin on each site. This is the ultra- 
strong coupling (i.e. local) version of a Mott insu- 
lator. Note that neither of these two states exhibits 
excitonic order. This can be seen by directly com- 
puting (a^b) = in either state. In fact, the opera- 
tors a^b and b^a act to transform the two phases into 
one another, i.e. move an electron from the lower to 
upper orbital or vice versa. 

2.3 Hopping terms 

To investigate further, we must introduce hopping 
between adjacent cells. We will principally consider 
the simplest such term, 



H , = ^2t{ata j + b\b j + h.c.), 



(6) 



where (ij) indicates that the sum is over nearest 
neighbor pairs of sites. Different hopping integrals t a 
and tb between the two orbitals can also be easily in- 
cluded, but do not change the results significantly, so 
we will keep t a = tb = t for simplicity (see however, 
the discussion of particle- hole symmetry in Sec. |5| 
surrounding Eqs. 51-55|). In general, there are also 
hopping processes connecting a and b orbitals. Due 
to the symmetry of the orbitals in Fig. 1, these occur 
only for next-nearest neighbors, 



H' 



t ab sign[(xi - Xj)(yi - yj)] 



a\b^ + h.c. 



(7) 



where the double angular brackets denote a sum over 
next nearest neighbors i and j. Note that the hop- 
ping matrix elements in Eq. |7| are real and vary in 
sign. The sign variations reflect the symmetry differ- 
ences (under rotations) between the s and d orbitals. 
The reality of the coefficients is a matter of conven- 
tion, which we fix by choosing the orbital wavefunc- 
tions to be real. We will assume, as appropriate in 
this example, that t a b <C t, so that H" is small, and 
can therefore be treated perturbatively. 

It is sometimes an important perturbation, be- 
cause it reduces the symmetry of the Hamiltonian. 
In particular, all of the terms in Ho + H\ + H' con- 
serve the number of a and b particles separately. 
Neglecting the a—b hopping, therefore, the model 
has SU(2) x U{1) x U(l) continuous symmetries, 
corresponding to conservation of spin, and a and b 
charges. The perturbation H" reduces the continu- 
ous symmetries of the model down to SU(2) x U(l) 
corresponding to spin and total charge, which are 
required by the physics of the system. Eq. [j] actu- 
ally still respects a number of discrete symmetries, 
such as b — > —b simultaneously with a ir/2 rotation. 
These symmetries, which in fact comprise the point- 
group operations of the square lattice, can be viewed 
as a residual discrete subgroup of the original U(l) 
present in the absence of H". We will see in the next 
section that this gives rise to an Ising symmetry un- 
der which the excitonic order parameters transform. 

3 Effective Theory for the Un- 
doped System 

In the central region of Fig. 2, e.g. for (U + V)/2 < 
fi < (U + 3V)/2, all sites are doubly occupied in 
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the strong coupling limit. Nevertheless, for 2Eq < 
U — V, the low energy states are highly degenerate. 
Well to the left of the thick vertical line, each a and b 
orbital is singly occupied, so that there are effectively 
two spin- 1/2 degrees of freedom in each unit cell. 
In the infinite coupling limit these are completely 
free, but they will interact due to virtual hopping 
processes when H' (and H") is included. Far to the 
right of the vertical line, the unique low energy state 
consists of a doubly occupied b orbital in each unit 
cell, and hopping is unimportant. As the vertical 
line is approached from either side, virtual hopping 
processes can induce interactions involving all five 
low energy states. 

3.1 Bosonic t—J model 

In this subsection, we develop an effective model for 
the interesting region near the vertical line. In this 
region, it is necessary and sufficient to truncate the 
Hilbert space to just the five low-energy states in 
each unit cell (although higher energy states must be 
kept in virtual processes). The physics is amusingly 
similar to a sort of generalized bosonic t—J model. 
On the left-hand side of the thick vertical line, each 
unit cell is occupied by two spins. At second order in 
H', these interact via effective exchange interactions, 

-^eff = ^ ] J ( ' Sja + Sib " Sjb J — ^ JllSia ■ Sift, 



(ij) 



(8) 



where J = t 2 /(V + 2E G ). This exchange constant 
may be obtained by computing perturbatively the 
energy difference between singlet and triplet states 
on a bond to second order in the hoppings, and ne- 
glecting the deviation from the vertical line (i.e. set- 
ting U = 2E G + V) in the denominators. The latter 
approximation is valid provided \ U — V — 2E G \ <C V. 
Well to the left of the vertical line (in particular 
when U — V — 2E G 3> t), no doubly-occupied b states 
are present, and Eq. || is a complete model. It de- 
scribes two ferromagnetically bulk coupled Heisen- 
berg spin-1/2 antiferromagnets. On a hypercubic 
lattice (square or cubic in two or three dimensions, 
respectively), one expects long-range antiferromag- 
netic order of spins on the same orbital sublattice, 
with a and b spins aligned parallel at each site. 

As the vertical line is approached, the energy cost 
of a doubly-occupied b orbital is reduced towards 
zero, and they must be introduced into the lattice. 



Unit cells with both electrons in the b orbital act as 
"holes", having no associated local moment. Unlike 
the usual t-J model holes, they are, however, bosonic 
and neutral (relative to the magnetic state, they rep- 
resent the removal of an a electron and replacement 
with a b electron). Hole hopping occurs at second 
order in t, : 



H. 



» (ij) 
+ Y. V hhh\W)h 3 (9) 

(ij) 



where Hh = 2E G -U + V + t 2 /(2V)-t 2 /[2(2E G + V)] 
is the hole "chemical potential", t^ = t 2 /(2V), 
V hh = t 2 /V - \t 2 /{2E G + V), and Vy = (§ + 2S ia ■ 
Sj a )(^ + 2Sib- Sjb) is the operator which interchanges 
the spin states at sites i and j. Like in a conven- 
tional doped anti-ferromagnet, the presence of the 
Vij operator in the hopping term leads to difficulties 
of hole motion in an antiferromagnetic spin back- 
ground. Naive successive hopping of a single hole in 
an antiferromagnetic state results in a generalization 
of the well-known "string" of misaligned spins in its 
wake. 

Introducing the a — b hopping term (Eq. [/]) affects 
the system in several ways. There are renormaliza- 
tions of the coupling constants in Eq. |9| and Eq. [§[ 
of order t 2 b /V, t 2 b /(V + 4Eg). Since, by assump- 
tion, tab «C t, these are negligible. New exchange 
couplings are also generated between next-nearest- 
neighbor a and b spins, which were not previously 
present. Because they are small, unfrustrating, and 
break no additional symmetries, these are also neg- 
ligible. The most important effect is to introduce a 
term which violates /i-particle conservation: 



Tjnnn 



m) 



a;ij 



Ofi^+h.c. 



(10) 



Here O^,^ creates a triplet of spin one states of 

a/b particles on the pair of sites ij, O^,^. creates a 
singlet of a/b particles on this pair, and the "fugac- 
ity" y = 2t 2 ab /V. Note that although Eq. [lC] violates 
conservation of the number of "holes" , it creates and 
annihilates them only in pairs. There thus remains a 
conserved Ising charge or parity (= ^ /ij/i i (mod2)), 
signifying whether the number of holes is even or 
odd. This parity can be traced back to fact that the 
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two orbitals on each site transform differently under 
spatial reflections. 

3.2 Pseudo-spin description 

To understand the behavior of this model, we now 
introduce a useful reformulation. Formally, the 
five possible states on each site can be viewed as 
different quantized values of a generalized pseudo- 
spin, and the above terms then take the form of 
nearest-neighbor interactions between these spins. 
In particular, we define five states per site via 
ID = a\bt\v),\2) = atb]\v),\3) = a\b\\v),\A) = 



Lastly, we consider the effects of the hole-pair cre- 
ation and annihilation terms in Eq. |l0|, which can 
also be transcribed into the pseudo-spin language. 



One finds H^g 



fl£ + 3$f, where 



Heft 



i 3 ' i j 



-35*7-35 



_^Tl5^7-45 _ ^7-45^7-15 1 q-b2q-b2 1 ^53-7-53 

1^ lj 1^ lj ~~v 1{ lj \ 1{ 1 j 



srtiLq'tyl _ /T-54/T-51 

' 3 i 3 



(12) 



1 5) = The Hamiltonian can be 

rewritten in terms of 5 x 5 spin matrices T^ u , where 
(fi r \T^\u') = 5^ 1 ' 5 UU ' . Neglecting for the moment 



the hole non-conserving terms in Eq. 1C, H£g 



H^ff + const., where 



TTPS. 

-"cff" 



J±< 



(ij) 



j}+ J z T*Tf-HY^ 

i 

(11) 



and2? = (£* 



Mi 



2tfii Jz 



The generalized ex- 
= Vhh, and Zeeman 



^4 

change constants J± 
field H = dV hh /2 - fi h . 

This form of the Hamiltonian exposes a strong 
similarity to the spin- 1/2 XXZ model in a Zeeman 
field. In particular, the "boson hopping" J± is 
analogous to an antiferromagnetic in-plane exchange 
{SfSj terms ), spin-boson interaction J z to an anti- 
ferromagnetic Ising exchange, and 7i to field. 
For J± 3> J z and TL not too large, one expects the 
analog of canted XY antiferromagnetism, while for 
Jz 3> Jx-, one expects instead z-axis Ising antiferro- 
magnetism up to a threshold value of \TL\. For large 
fields, \TL\ 3> J±,J Z , one expects ultimately fully 
polarized states, which correspond to the Mott and 
band insulators for TL > and TL < 0, respectively. 

Surprisingly, H^g displays an enormous SU(4) 
invariance under T 5 ^ -► *£t=i u nuT hv , -► 

Ej=i U^T vS , where u is an su ( 4 ) matrix - su ( 4 ) 
symmetry is expected to be a good approximation 
over a range of energies, because in the physical limit 
V < U ~ E G , J H < J < J±,J Z ,H. Thus we will 
take the approach of first solving the SU(4) invariant 
model, and considering successively the exchanges 
J and Jh, which reduce the symmetry of H c g to 
SU(2)xSU(2) (independent physical spin rotations 
of the a and b moments) and SU(2)xU(l), respec- 
tively. 



The coupling Jj oc y. While it is perhaps not com- 
pletely transparent in this notation (a better nota- 
tion for this term will be introduced in next subsec- 
tion - see Eq. |l9|) , the effect of is to further break 
the SU(2)xU(l) symmetry down to SU(2)xZ 2 . The 
Z2 invariance is the remnant of the physical parity 
symmetry discussed in the previous subsection. 

T z 

3.3 Mean-field theory and undoped 
phase diagram 

We expect that a simple Weiss mean field theory 
(MFT) gives qualitatively correct results for the sto- 
ichiometric phase diagram, as it does for the ordi- 



nary XXZ+Zeeman model. [16| Neglecting H^ s , the 
MFT consists in replacing 



\ i / i ' i \ j 



<£t5 Irp5fl.\ 



(13) 



for each bond i,j on the lattice, and similarly for the 
TfT? interaction. With this replacement, the Hamil- 
tonian decouples on different lattice sites, and the 
problem reduces to solving self-consistently the ap- 
propriate single-site problems. As an antiferromag- 
netic solution is expected, this amounts to equations 
for the (8 component) transverse staggered magneti- 
zation, defined by (T/* 5 ) = (— l) l [n 2 _^~ + in 2 ^], and 
the uniform and staggered z-axis magnetizations, de- 
fined by (T?) = m z + (-l)W Because of SU(4) 
symmetry, all orientations of are degenerate, and 
it is sufficient to assume = n±5 kl . In this sub- 
space, the equations of MFT become identical to 
those of the conventional spin-1/2 XXZ antiferro- 
magnet in a Zeeman field. These equations were 
solved in Ref. [17|. The resulting phase diagram is 
shown in Fig. |. 

Since J± > J z , we expect transverse pseudo- 
spin polarization, (T^ t5 ) 7^ 0, provided \TL\ < TL C = 
dj±. Remarkably, the transverse components of the 
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Figure 4: Mean-field phase diagram of the fully 
SU(4)-invariant pseudo-spin model. A large pseudo- 
Zeeman field (which scales linearly with the orbital 
splitting Eq) stabilizes either the Mott insulating 
(MI) or band insulating (BI) state, depending upon 
its sign. For "in-plane" anisotropy (J7j_ > J z ), the 
intervening phase is an excitonic insulator (EI). In 
the opposite limit ( "Ising" anisotropy) it consists of 
a micro-phase-separated state with a checkerboard 
pattern of alternating band and Mott insulating con- 
figurations at the lattice scale. In obtaining the 
pseudo-spin model from the strong-coupling limit of 
Eqs. £||], one finds in-plane anisotropy, and the inter- 
mediate state is expected to be excitonically ordered. 

pseudo-spin operator are exactly the excitonic or- 
der parameters. In particular, straightforward alge- 
bra shows T^ 5 = (— A|j_, Ajf, — Aj_j_, Aj_j). Thus for 
J±_ > J z , MFT predicts an excitonic insulator. 

We now turn to the evolution of the ground state 
in this regime on introducing the symmetry-breaking 
terms in H^. In their absence, the excitonic or- 
der parameter can "point" in any direction which 
is equivalent under the broken SU(4) symmetry. 
Within MFT, this amounts to complete freedom to 
choose the four complex components of A a p, sub- 
ject to the constraint Tr A 1 " A = \(l - H 2 /H 2 C ) = 
Aq/2. In term of singlet and triplet components 
defined by Eq. 0, this constraint simply implies 



A s r + A? -A t 



A 2 



The perturbations in H*L can 



be viewed as "anisotropics" favoring sub-manifolds 
within this space. 

To clarify the nature of the anisotropy terms, it 



is helpful to work with the mean-field wavefunction, 
l*o> = Ui£}\BI), where \BI) = Ui^v) is the 
non-interacting band-insulating state, and 



e\ = c\ i + (-i) 1 



af3 



(14) 



is a local "exciton creation operator". Here |c| 2 = 
(1 — 7i/7i c )/2. It is now straightforward to evalu- 
ate the expectation value of H£ s in the mean-field 
ground state. Up to a constant for fixed TrA^A, 
on a hyper-cubic lattice one finds the bulk energy 
density 



e b = L- d {H s cS ) = 2JTr(A t A) + J H |Tr A 



(15) 



where J = da- d J/2\c\ 4 and J H = aT d J H /2\c\ 2 . The 
above terms are essentially completely determined 
by the SU(2)xSU(2) and SU(2)xU(l) symmetries. 
To proceed, we employ two identities derivable from 
Eq.g 



Tr 



(AtAY 



Tr A 



1 

8 
1 

8 

A,. 



a:a s + A* • A, 



A* s A t + A s A* t -iA t AA* t (m 

(17) 



By assumption, J Jh, so that the first term 
in Eq. [l5| creates the dominant splitting of the 
SU(4) ground-state degeneracy. The low-energy 
sub-manifold thus consists of the order parameters 
which minimize Tr(AtA) 2 . Eq. ^ then implies 
A* s A t + A S A* - iA t A A* t = (note that the first 



term in Eq. 16 is constant and equal to Aq/8). The 
physical content of this condition is made clear by 
calculating the mean spin polarization on the a site 



using the mean-field wavefunction in Eq. 14 

s a = (S a ) 



v iA t A A* + A* A t + A s Aj? 



(18) 



Thus the influence of the exchange coupling J is to 
favor states with s a = 0. 

This condition still allows a fairly large range of 
states, the simplest of which are pure singlet (|A S | = 
A , A t = 0) and pure collinear triplet (A s = 0, 
At 7^ 0, A t A Af = 0) orderings. The additional 
effect of the Hunds-rule ferromagnetic coupling Jh 
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is to introduce a small extra "mass" for the singlet 
order parameter, favoring a pure triplet state. 

The phase of the triplet order parameter is deter- 
mined by the "Ising anisotropy" terms in Eq. To 
see this, we rewrite T m5 and T 5 ^ directly in terms of 
A. One finds 



H. 



eff 



J/Tr A, A 



+ a\a] 



(19) 



Note that Eq. 16 explicitly breaks the U(l) symme- 
try of phase rotations of A, down to the Ising invari- 
ance A — > — A. If Hi* is considered a weak perturba- 
tion, it can be treated by simply evaluating its expec- 
tation value in the mean- field ground state (Eq. [14]), 



giving (Ai 



(A? 



: A, since i and j are next- 
nearest neighbors. Using Tr A 2 = (A 2 + A t • A t )/2, 
one finds (since J 1 > 0) that Eq. 19 favors an imagi- 
nary triplet order parameter At = —A*. This is dif- 
ferent from the weak-coupling treatment of Ref . Q , 
in which a reaJ triplet order parameter was found to 
be preferred. 

Unlike in superconductivity, the phase of the ex- 
citonic order parameter has physical significance, as 
discussed by Halperin and Rice[||]. In particular, it is 
straightforward to show that a real At order param- 
eter corresponds to a non-zero average spin density 
within the unit cell of the crystal, while for At imag- 
inary, the spin density is zero but there are instead 
non-zero spin currents. The imaginary triplet state 
obtained here is therefore a sort of spin "flux phase" 
with non-zero spin currents. See Sec. || for a more 
in-depth discussion. 

Apart from this difference, the strong-coupling re- 
sults of this section are in very close agreement with 
the weak-coupling results of Ref. ||. Indeed, not too 
much significance should be attached to the differ- 
ence in phase of the order parameters, as indeed the 
models are in any case not completely identical. In 
fact, the detailed correspondence of results up to this 
point strongly argues for a continuous smooth inter- 
polation ("adiabatic continuity") of most physical 
properties of such systems as the overall interaction 
strength is increased from small to large values. 

Finally, we comment on the modifications to the 
SU (4) -invariant phase diagram in the presence of the 
symmetry-breaking terms in Eqs. ||, [l(]. As argued 
above, these favor an imaginary triplet state when 
TL = 0. Inside the Mott insulator, these terms stabi- 
lize an antiferromagnetically ordered magnetic state. 
On approaching the Mott insulator boundary, there- 



fore, we expect the emergence of magnetic ordering. 
This implies the existence of at least one additional 
phase boundary separating the triplet EI (which has 
no non-zero spin density) from a magnetically or- 
dered EI with non-zero average spin density, some- 
where inside the region in which the EI phase occurs 
in the SU(4)-invariant model. 



4 Doping 

In this section, we consider the behavior as a low 
density of electrons is added to the system. In the 
strong-coupling limit, this reduces to an effective 
i-J-like model, in which the Hilbert space is re- 
stricted to states in which all sites (unit cells) are ei- 
ther doubly (corresponding to the excitonic pseudo- 
spins modeled above) or triply occupied, the lat- 
ter containing one a and two b electrons. The sys- 
tem is then governed by an effective Hamiltonian 
Hdope = H^g + + VH'V, where V projects onto 
this restricted Hilbert space. 

As many years of work on high-T c superconduc- 
tivity has taught us, the problem of doping a cor- 
related (Mott) insulator, particularly with spin (and 
here pseudospin) ordering, is extremely complex and 
difficult. Here, we will adopt the absolute simplest 
approach extending the above MFT to the low elec- 
tron density limit. We assume, as suggested by the 
weak-coupling analysis, that the essential ingredient 
for excitonic ferromagnetism is the approximate en- 
hanced (in this case SU(4)) symmetry of the effective 
Hamiltonian. In considering the doped state, then, 
it is crucial to determine in what way the added 
electrons affect the splitting of the degenerate SU(4) 
ground-state manifold. 



4.1 Variational treatment for a single 
electron 

In the strong-coupling limit, the majority of the 
energy of an added electron is kinetic, since t 3> 
J ', J_\_ ~ t 2 /V, etc.. Just as in the simpler but much 
studied t-J model for the cuprates, coherent motion 
of an added electron, however, is greatly hindered 
by (pseudo)-spin ordering of the insulating back- 
ground. Moreover, coherent motion is possible to 
a varying degree depending upon the precise nature 
of the background. We first consider this effect for a 
single added electron using the variational method. 
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A natural variational ansatz is 

id j^i 



(20) 



where both the doped electron's wavefunction ipi a 
and the excitonic order parameter A a p (implicit in 
£j) are considered as variational parameters. For 
fixed Tr At A, the energy depends only upon J, Jh, 
and t. In particular, one finds 

e v = L-^il^dopel*!) = e b (l - 2d(a /L) d ) + L~ d e e 

(21) 

where ao is the lattice spacing, 

ee = tJ2^aTa^ j/3 , (22) 

and the matrix T a p = \c\ 2 5 ol p+\c\~ 2 (AA^) . Phys- 
ically, we identify the first term in Eq. |2l] as the 
bulk energy density, reduced by the presence of a 
single doped electron (occupying the volume frac- 
tion (a/L) d ). In the second term, the quantity e e is 
then readily interpreted as the energy of the added 



electron. Eq. 22 is then a hopping Hamiltonian for 
this electron. In a polarized excitonic background, 
this hopping is in general non-diagonal in spin. In 
terms of singlet and triplet components, 



:i + \c\- 



+ s a ■ a 



(23) 



where s a , the mean spin polarization on the a site, is 



given by Eq. 18. Minimizing Eq. [2^ in the space of 
normalized wavefunctions ipi a gives the tight-binding 
Schrodinger equation, 



(ji) P 



(24) 



where the angular brackets indicate a sum over the 
nearest neighbors j of site i. The single-particle 
eigenstates of this equation are plane waves with 
spins polarized parallel and antiparallel to s a , with 
eigenvalues 



e e ±(k) = 2t 



1 + \c\ : 



± 141 



[cos kia ] , (25) 



i=i 



where ao is the lattice spacing. The location of the 
minimum-energy electronic excitations depends cru- 
cially on the magnitude of s a , and hence 7i. When 



H > Tic/2), electrons with spin parallel and antipar- 
allel to s a have minimum energy at different points 
in momentum space. Such large values of H cor- 
respond to strongly overlapping bands, close to the 
boundary between the Mott and excitonic insulators. 
For simplicity, we will specialize to the case when 
\s a \ < (1 + |c| 2 )/2, which occurs for H < H c /3. In 
this case, the minimal energy single-particle energy 
excitations for both spin orientations have momen- 
tum k = (vr, • • • ,7r). Furthermore, the optimal spin 
orientation is parallel to s a . Such an electron takes 
advantage of the "Zeeman" energy due to the ex- 
change field (proportional to s a ) generated by the 
"core" spins (i.e. the spins of the two electrons per 
unit cell present in the insulator). 

In the undoped system, however, s a = 0, due to 



the anisotropy in Eq. 15. We therefore expect that 
the optimal order parameter in the doped system 
is determined by a competition between these two 
terms. With some algebra, it is straightforward to 
verify that, due to the Hunds-rule term Jh, the com- 
plex pure triplet state (i.e. with A ( A / but 
A s = 0) is always more energetically favorable than 
a singlet-triplet coexistence (with Re A*A t ^ 0).[f[g|] 
Without loss of generality, it is thus convenient to 
choose a spin quantization axis, letting 



At = Aq(cos 6>x + i sin 6y), 



(26) 



One then finds s a = — (Aq/2|c| 2 ) sin(2#)z. In any 
such state, Sb = s a , so that the core spins also con- 
tribute to the ferromagnetic moment. 

4.2 Free Fermi gas approximation 

It remains to determine the optimal angle 8. To 



proceed, we need to extend Eq. 21 to a small but 
non-zero density of doped electrons. At low den- 
sities, it seems natural to neglect interactions be- 
tween doped electrons, and use the simplest possi- 
ble free Fermi gas estimate for the electronic dopant 
energy. In particular, we approximate the energy of 
the system as the sum of two contributions: a "bulk" 
contribution from the undoped unit cells containing 
two electrons and a spatially uniform order parame- 
ter A a p, and a "dopant" contribution, approximated 
by the energy of a free Fermi gas of electrons with 
dispersion given by Eq. |25|. For concreteness, the 
detailed formulae are presented in the following for 
three spatial dimensions (d = 3). At low densities, 
only single-particle states near k = tv = (tt,tt, tt) are 
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occupied, so it is convenient to expand around this 
point, k = 7r + q, yielding the dispersion 



e e ±(q) 



-2t 



(l + |c| : 



± 141 



(27) 



Here we have re-instated a (shifted) chemical poten- 
tial jx to control the density of doped electrons. It 
is both convenient and physically helpful to work 
at fixed chemical potential rather than fixed charge 
density, as this allows naturally for the possibility 
of phase separation. As is perhaps not surprising 
based on the results of weak-coupling analysis || JtJ, 
we will see that phase separation does indeed occur 
in a physically interesting parameter range of the 
model (at least within this approximation). 

Because we are interested in the energy density 
only insofar as to determine the angle 9, we neglect 
in the following all terms independent of 9. Inserting 



Eq. 2£ into Eq. 15 gives the bulk energy 



e b = [3JA^/(8ag|c| 4 )] sin 2 29 + const. 



(28) 



(in three dimensions). This must be added to the 
ground state energy of the free Fermi gas of doped 
electrons. Simple but tedious algebraic calculations 
lead to the final expression for the total energy den- 
sity of the system: 



I5 2 



-£(<?, 7) 



(29) 



where 



2\2 



3J(l + |c| 2 ) 
8a 3 
A 2 

|c| 2 (l + |c| 2 )' 
5tt 2 (1 + |c| 2 ) J 
16\/6 t 



9 
A 

7 



sin 29 
1 + 



3t(l + |c| 2 ) ' 



X/8. 



(30) 
(31) 

(32) 
(33) 
(34) 
(35) 



The function £(17,7) is straightforwardly related to 
the energy density of the three-dimensional free elec- 
tron gas in a Zeeman field. In general it depends 
not only on g and 7, but also on 5. For simplicity, 
we will assume \5g\ <C 1, which holds near to the 



excitonic insulator-band insulator boundary, and is 
satisfied more generally in the interesting region of 
the phase diagram (where 5 is 0(S C ), since S c <C 1 
in the strong coupling limit J/t <C 1 — see Fig. ||). 
In this case, £ (5,7) becomes independent of 5. Its 
functional form is 



£(9,1) = ]T(7 + ^) 5/2 e( 7 + ^), 



(36) 



=±1 



where 0(7) is the Heavyside step function. Eqs. 29- 
pq give the energy density of the system func- 
tion of chemical potential jl (through 7) and order 
parameter angle 9 (through g). The optimal exci- 
tonic angle 9 is determined by minimizing e(/2, 9) at 
fixed jl. If jl and 9 are known, the density of doped 
electrons x and itinerant magnetization density ma 
are then given by the free-fermion results: 



V2 
6n 2 a 3 



Y^{l + zgf/ 2 @{ 1 + zg), (37) 



z=±l 



^y>( 7 + 2a) 3/*e (7 + zrf 38) 



Note that the system is doped (i.e. i/O) whenever 
7 + \g\ > 0. One should keep in mind also that the 
full magnetization density m = mit+m core includes a 
contribution m core = — (Ag/|c| 2 a 3 ) sin(2#) from the 
core spins. 

Eqs. 29-3^ completely determine the state of the 
system at zero-temperature as a function of A and 
5. The mathematical problem of minimizing ej is al- 
gebraically quite tedious, and significant care must 
be taken to avoid spurious local minima and saddle 
points. The results of a careful study are shown in 
Fig. H All the phases shown are excitonically or- 
dered, but differ in doping x, excitonic angle 9, and 
magnetization m. The properties of each are sum- 
marized in Table [j]. In the strong-coupling limit, 
we expect (see Sec. ||) 5/5 c 3> 1, in which case 
there is a direct Erst-order transition from a param- 
agnetic excitonic insulator (EI) to a fully-polarized 
ferromagnetic metal (FMFP*). The phase bound- 
aries in Fig. H variously indicate first (discontinuous) 
and second(discontinuous) order transitions. All the 
vertical phase boundaries denote continuous transi- 
tions, while most of the transitions on curved phase 
boundaries are discontinuous. The exceptions are 
the PPFM*-FPFM* boundary (which is everywhere 
second order) and the lower-portion of the FPFM- 
FPFM* transition line, which is continuous below 
the tricritical point indicated in the figure. 
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Figure 5: Free- Fermi gas approximation to the zero- 
temperature phase diagram. The abscissa X/5 C is 
the non-dimensionalized chemical potential, and the 
ordinate S/5 C indicates the strength of excitonic or- 
dering relative to super-exchange interactions (see 
Eqs. ||-|| for the precise definitions). The proper- 
ties of the various phases shown are listed in Table |l[ 
For larger A (not shown), both excitonic ordering 
and ferromagnetism disappear, owing to reduction 
of the amplitude of A. 

Which portion of this phase diagram is most 
physically significant? In the strong-coupling limit, 
5 C <C 1, and it therefore seems reasonable to suppose 
6/6 c 3> 1, so that the system undergoes a simple and 
direct first order transition from the undoped and 
paramagnetic EI to the fully-rotated half-metallic 
ferromagnet, FPFM*. Coincident with this transi- 
tion is a jump in the electronic charge density x, 
from zero in the insulator to a non-zero value in the 
metal. 



phase 


doped? 


mag. 


pol. 


angle 


EI 


no 


para 


n/a 


6 = 


PM 


yes 


para 


n/a 


6 = 


PFM 


yes 


ferro 


partial 


< 6 < vr/4 


FM 


yes 


ferro 


full 


o < e < vr/4 


PFM* 


yes 


ferro 


partial 


6 = vr/4 


FM* 


yes 


ferro 


full 


6 = vr/4 



Table 1: Phases of the doped EI in the free- Fermi 
gas approximation. The five columns list the ab- 
breviation, presence or absence of doping, magnetic 
order, degree of polarization, and excitonic angle (6), 
respectively, for the six phases. 

greater than two per unit cell. Nevertheless, if, as 
supposed, SU(4) symmetry is a good approximation, 
then many other possible states must necessarily be 
nearly as low in energy. In the hope that the truth 
may ultimately be decided by experimental measure- 
ments, it seems useful to delineate the physical char- 
acteristics of each of these phases. 

With the exception of the non-collinearly ordered 
states, the analysis of the next few paragraphs is 
identical (though in somewhat different notation) to 
that of Halperin and Rice||. First, let us consider 
the existence of a time-averaged magnetic moment. 
In the tight-binding formulation, the electron field 
operator is expanded in terms of Wannier orbitals, 



en. 



+ <j) b (r - Ri)6j, 



(39) 



where <p a / b (r) is the Wannier function for the a /b or- 
bital, and we neglect the other (unoccupied) states. 
Consider next the spin density operator. We will 
assume for simplicity (though this is not essential) 
that each Wannier function has support only within 
one unit cell. Eq. |3^ then leads to a representation 
for the spin density operator S, 



5 Discussion 

5.1 Symmetries and properties of exci- 
tonic insulators 

The model introduced in Sec. || contains many pos- 
sible excitonically ordered states in various regions 
of its phase diagram. In the undoped case, we 
have argued that a simple paramagnetic collinear 
triplet ordering is most likely, while a state with 
At A A* ^ obtains for electron densities slightly 



2S(r) = ip^aip, (40) 
= \Mr)\ 2 (^aa) + \Mr)\ 2 (tfab) 
+ 0*(r)&(r)(at<?&) + #(r)^(r)(&t;? ( flLL) 

To proceed, we choose both Wannier functions to be 
real. Then for the undoped case, the spin density 
can be rewritten in terms of s a i b and Aj: 



Six) = \<t> a {v)\ 2 s a + \<t> b (r)\ z s b + 20 a 



2-> 



6 b (r)ReAi. 

(42) 
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Recall further Eq. O and its analog for s&: 



rate, Eqs. nU2| then imply that 



?a/b 



4c 



iA t A A* ± 



A: A t + A S A* 



(43) 



There are thus non-zero static local moments when- 
ever Re A* A<, ImAf A A*, or Re At are non-zero. 
In the simplest such states, A t = |A t |e, where e is 
a real unit vector. In this case, there is a spatially- 
varying static moment within the unit cell oriented 
along the e axis. The net moment (integrated over 
the unit cell) is, however, zero, unless A t A A£ 7^ 0, 
in which case the real and imaginary parts of At 
are both non-zero and not parallel. In addition to 
the net ferromagnetic polarization along ImA^A A?, 
such states have a non-collinear static spin density 
in the unit cell. The net moment along these other 
directions remains zero. To see why such states sus- 
tain a net polarization, consider the particular case 
given in Eq. ^6|, with excitonic angle 6. One can then 
use Eq. [2] to rewrite the order parameter matrix as 



A = -y [(cos 6 + sin d)a + + (cos 9 - sin 9)a~ 



(44) 



Inspection of the mean-field wavefunction, Eq. 14 



and Eq. 44 immediately shows that the amplitude 



for up and down spins are unequal, so long as 6 is 
not a multiple of ir. 

Some confusion may arise in the reader with re- 
gard to time-reversal symmetry. It appears surpris- 
ing to have At and iAt A Ajf, the latter containing 
a cross-product, both contributing to s a /b- In fact, 
both terms transform like a spin under time-reversal. 
This is simplest to see in the path-integral represen- 
tation of the quantum system, in which the Fermion 
operators are replaced by time-dependent Grassman 
fields a a — > a a (t), a a — * a a (t), and similarly for 
b a ,b a . The Grassman fields then transform under 
time-reversal according to 



a a (t) 
a a (t) 

b a (t) 
b a (t) 



-<,&/*(-*)■ 



(45) 
(46) 
(47) 
(48) 



Note the important minus sign in the above trans- 
formation, which is possible because a a , a a (b a , b a ) 
are independent fields (not related by complex con- 
jugation) in the path integral. This reflects the anti- 
unitary nature of time-reversal symmetry. At any 



A s 
A f 



-A* 



(49) 
(50) 



under time-reversal. The combination of complex 
conjugation and the minus sign for At imply that 



both terms in Eq. |43|are odd under time-reversal. In- 
deed, the necessary and sufficient conditions for bro- 
ken time-reversal symmetry is Re (Af) 7^ and/or 
Im(A s ) ^0. 

A perhaps surprising consequence of Eqs. 41 is 
that apparently if ReA*A t / but Im A t A A| = 0, 
there is no net magnetization. In fact, this result ap- 
plies only to the particular undoped model consider 
here, and is a consequence of a special variety of 
particle- hole symmetry (which we denote PH). To 
make this explicit, define a hole creation operator 



Then the electron number operator can 



be rewritten as 



n 



a)a + b [ b = 2 + a)a-b%. 



(51) 



In the undoped system, the mean number of elec- 
trons per unit cell is two, so that (a^a — ¥b) = 0. 
Thus precisely at this density, and only at this den- 
sity, we may entertain the possibility of symmetry 
under the transformation PH: 



b a ~^ph a a . 



(52) 



In the new variables, the excitonic order parame- 
ter becomes A = a} a a v a ^^. Thus A — > — A T (here 
the superscript T indicates the matrix transpose) 
under PH. Also useful is the operator a)a — tfb = 
a) b + ¥b — 2 (proportional to T z in the n = 2 sub- 
space), which is invariant under PH. Thus H^g (see 
Eq. [IT]) is PH- invariant. Similarly, it is straightfor- 
ward to show that under PH, the two spin operators 
are exchanged: 



Sa 



Sb- 



(53) 
as can be 

easily shown. Thus the undoped Hamiltonian is in- 
variant under PH. Considering the order parameters, 
we find that 



Thus is also PH- invariant 



as is H* s , 



At 



>PH 



A, 



-At 



(54) 
(55) 



Thus the combination Re A*Af is odd under PH, 
and hence cannot give rise to a total moment, since 
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Stot is PH-invariant. Im At A A*, however, is PH- 
invariant, and can hence couple directly to a ferro- 
magnetic moment. 

It should be stressed that PH is not a microscopi- 
cally exact symmetry, even in the stoichiometric sit- 
uation. It occurred in the above analysis only be- 
cause of the arbitrary choice of equal hopping be- 
tween a and b orbitals, t a = % = t, in Eq. ^. In 
general, one expects t a / tb, which leads to differ- 
ent anti-ferromagnetic exchange constants between 
a and b spins in H^ s , Eq. ||. Different exchange con- 
stants destroy the invariance of the Hamiltonian un- 
der the interchange of a and b spins, Eq. [5^, which is 
the effect of PH. It is straightforward to show that, 
when this asymmetry is included in the microscopic 
Hamiltonian, states with Re A* At 7^ are also fer- 
romagnetic. In addition, even in the model with 
t a = tb, doping breaks the PH symmetry, and gives 
rise to a ferromagnetic moment in the Re A* At 7^ 
state. 



Considerations similar to those above Eq. 41 apply 



to the electronic charge density (/?), current density 
(I) , and spin current density (J^ v ) operators. One 
finds 



p(r) 
f(r) 



-e\4> a {r)\ 2 n a - e\(p b (r)\ 2 n b 
-2e0 a (r)0 6 (r)Re A s , 



-Im A, 



m 



a (r)V<&(r) - b (r)V^ a (r) 

>(r)e^ a (r)fl5.6) 
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2m 



In the final equation above, J^ u is the current den- 
sity for spin polarized along the p axis propagating 
in the v direction. For completeness, the mean-field 
expressions for the number of a and b particles are 



(a+a) 



= |c|- 2 TV A f A, 
2 - n a . 



(57) 
(58) 



From Eqs. 56, we can read off the physical inter- 
pretation of the various other types of ordering. If 
Im At 7^ 0, there is a spontaneous spin current in the 
unit cell. This is necessarily the case for any state 
with At A Ajf 7^ 0, which, as discussed above, also 
exhibits non-collinear static moments. The simpler 
state with At = i|At|e has only the spin currents, 
and is the magnetic analog of a "flux phase" in mod- 
ern terminology. Similarly, if the singlet order pa- 
rameter has an imaginary part ImA s 7^ 0, there are 
non-zero charge currents within the unit cell. This is 



exactly a flux phase. Finally, a real singlet order pa- 
rameter, As = A* 7^ 0, gives rise to a charge-density 
p(r) that breaks the point group symmetry of the 
crystal, since (j) a (r)(J)b(r) is not a scalar. 

Another importance characteristic of the phases 
with triplet ordering is a finite (transverse) uniform 
spin susceptibility. This is a very general conse- 
quence of broken spin-rotational invariance. In the 
simplest collinear triplet states, At = Aoe*^e, where 
e is a real vector. The elementary excitations of the 
symmetry-broken state can then be classified only 
by their spin along the triplet axis, Stot • e. The 
transverse components of Stot> however, do not 
commute with At. An applied Zeeman field along 
one of these axes therefore immediately acts to mix 
together the former ground and excited states. It is 
fairly straightforward to demonstrate by this mecha- 
nism a constant transverse spin susceptibility for the 
collinear ly-order triplet states. For non-collinear ly 
ordered triplets, we conjecture that all components 
of the uniform susceptibility are finite. This distinc- 
tion is most likely primarily academic, as experimen- 
tally available samples would presumably break up 
into domains with random orientations of At, thus 
effectively isotropizing the bulk susceptibility. Very 
crude estimates for the magnitude of x can be obtain 
in both the strong and weak coupling limits of exci- 
tonically ordered states. In strong coupling, the sus- 
ceptibility can be computed by naive perturbation 
theory in the mean-field ground state. In the opti- 
mal case (TC = 0), in which the excitonic ordering is 
maximal, one finds x ~ I^b/Jl-, where J± ~ t 2 /V 
is the characteristic stiffness for excitonic ordering 
(see Sec. |||), and p B is the Bohr magneton. In weak 
coupling, the susceptibility is approximately equal to 
the free-electron value, x ~ D(ep)p B , where D(ep) 
is the density of states at the Fermi energy. 

A puzzling aspect of the experimental data on the 
hexaborides is the absence of a substantial gap in 
optical conductivity measurements in the undoped 
materials. In general, the excitonically ordered in- 
sulators discussed here are expected to exhibit hard 
optical gaps (i.e. complete absence of weight in cr(u>) 
at small u) at low frequencies and zero temperature, 
so this is an important point which such a theory 
must contend with. Several possible physical situ- 
ations can, however, resolve this apparent discrep- 
ancy. In the weak and intermediate coupling limits, 
it is possible to sustain a metallic state simultane- 
ously with excitonic order. This requires imperfectly 
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nested Fermi surfaces - a detailed investigation of 
this possibility is underway. Even in the strong cou- 
pling limit, it is also possible that a gap exists but 
is anomalously small. Indeed, in the present model, 
the optical gap can be estimated by considering the 
energy cost required to transfer an electron from one 
unit cell to its neighbor. At the optimal conditions 
for excitonic ordering, one has U = 2E G + V and 
a straightforward calculation of energies from Eq. ||] 
gives the optical gap A G sa V. Thus one expects at 
least A D <C U, E G . Note that in the strong-coupling 
limit, there is no universal relation between A G and 
the excitonic order parameter. 

5.2 Relevance to the hexaborides 

The models and discussions in this paper demon- 
strate the feasibility of a strong-coupling approach 
to excitonic ordering. A direct application of the 
results to the hexaborides is, however, not appropri- 
ate, due to the simplified nature of the Hamiltonian 
discussed here. It is possible to generalize the tight- 
binding model discussed here to a "two-band" (p 
and d orbital) Hamiltonian which more accurately 
models the physics of these materials. This model 
contains the significant new ingredient of orbital de- 
generacy, and hence considerable additional richness. 
Such orbital degeneracy is the tight-binding ana- 
log of the valley degeneracy encountered in band- 
theoretic treatments. The methods of this paper, 
however, remain applicable in this case as well. A 
thorough treatment of this problem presents an at- 
tractive and challenging theoretical opportunity. 

It is reasonable to ask at this point whether there 
are any experimental consequences of the excitonic 
scenario which are relatively model- independent, 
and can therefore be firmly stated in advance of more 
accurate results? For this, we look to the discussion 
of the previous subsection, focusing particularly on 
the properties in the undoped material. The exci- 
tonic scenario postulates symmetry breaking even 
without doping, which distinguishes it from, e.g. 
low-density ferromagnetism ala Wigner. All calcu- 
lations so far appear to favor triplet ordering, which 
implies first a constant (temperature independent at 
lowT) susceptibility in the insulator. Second, triplet 
ordering necessarily gives rise to either static spin 
moments or static spin currents (or both) within the 
unit cell. Because the latter are presumably difficult 
to observe, this is a less strong condition. Third, 
because the triplet state breaks spin-rotational in- 



variance, it implies the existence of two low-energy 
collective "magnon" modes (presumably dispersing 
as uj as v s \k\ at low energies), which could be ob- 
servable via inelastic neutron or Raman scattering. 
Fourth, an excitonic explanation for ferromagnetism 
upon doping requires that the "pseudo-spin flop" 
phenomena occur, and hence (in this sense) approxi- 
mate SU(4) symmetry. This approximate symmetry 
implies the existence of additional collective modes 
with small excitation gaps. 

It is natural to ask, given the above emphasis on 
the undoped state, whether the excitonic ferromag- 
net is itself truly a distinct phase of matter separate 
from the more familiar (theoretically) Wigner ferro- 
magnet? The answer depends upon the extent to 
and manner in which the dopant ions influence the 
behavior of the electrons. In the models investigated 
to date, the dopants influence the material only inso- 
far as to donate extra charge carriers, providing no 
perturbation to the potential felt by the electrons 
except to slightly increase the neutralizing positive 
background charge. In this treatment, the lattice 
point group symmetries are strictly maintained, and 
the excitonic ferromagnet is indeed a distinct state 
of matter: it exhibits more broken (point group) 
symmetries than the Wigner ferromagnet. In real- 
ity, the dopant ions most likely distribute randomly 
throughout the crystal, and thereby perturb the po- 
tential experienced by the electrons. This random 
potential explicitly breaks the lattice invariances, 
and washes out this sharp distinction between the 
excitonic and Wigner ferromagnets. 

Whether this effect is of practical importance is 
unclear. The large increase of conductivity upon 
doping suggests that the electrons are not strongly 
scattered by the dopant ions. In any case, the 
physics at low electron density is quite subtle. In 
particular, the environment around a nearly isolated 
dopant atom retains a large fraction of the symme- 
tries of the pure lattice (e.g. a Lanthanum dopant re- 
placing strontium preserves cubic point group sym- 
metries around the lanthanum ion). Because elec- 
trons interact with only one impurity at a time at low 
densities, the symmetry of the local environment is 
expected to improve the distinction between Wigner 
and excitonic ferromagnets. 

Clearly, further predictions are possible within 
more specific models. Several authors g, [JJ have 
recently pointed out the likelihood of phase sepa- 
ration at low electron densities. This occurs natu- 
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rally in the pseudo-spin flop picture, but because it 
has already been discussed, we will not dwell on it 
here. Probably most importantly, any excitonically- 
ordered state by definition breaks the point-group 
symmetry of the lattice. This symmetry-breaking 
is directly observable, but unfortunately depends in 
detail on the way it occurs. In particular, many of 
the triplet states that appear to be favored have less 
obvious order parameters, so that more work needs 
to be done to ascertain the appropriate experimental 
probes. Further modeling using the strong coupling 
approach promises to help resolve these and other 



issues. 
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